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SECTION  M 


INTRODUCTION 

This  report  presents  final  results  from  a  FY-89  VLF/LF  Data  Analysis  task  initiated 
in  support  of  the  VLF/LF  propagation  modeling  program  being  developed  by  the  Naval 
Ocean  Systems  Center  (NOSC).  Specifically,  this  paper  gives  the  results  of  comparisons 
between  VLF  signal  measurements  recorded  during  nighttime  hours  and  NOSC  model 
predictions  for  similar  nighttime  propagation  paths.  Measured  data  was  recorded  in 
1985  and  1986  aboard  the  merchant  ship  G'^'S  Callaghan,  and  daytime  propagation  com¬ 
parisons  against  that  data  were  documented  in  the  NOSC  report  entitled  "Evaluation 
of  the  NAVOCEANSYSCEN  Long  Wavelength  Propagation  Capability  Using  VLF  Data 
Collected  Aboard  a  Ship”,  01  November  1987  (termed  herein  the  Callaghan  Report). 
The  present  report  is  offered  in  satisfaction  of  CDRL  A002,  under  D.  O.  0046,  Contract 
N66001-88-D-0033. 

For  many  years,  NOSC  has  developed  computer  programs  and  models  for  the  analysis 
of  VLF/LF  propagation,  and  has  applied  this  capability  to  a  number  of  applications  for 
the  Navy  and  for  other  DoD  agencies.  These  past  efforts  have  utilized  separate  pro¬ 
grams  for  various  parts  of  the  calculations  and  have  generally  been  performed  by 
experienced  NOSC  personnel  capable  of  judging  the  correctness  of  resulting  outputs. 

More  recently,  NOSC  has  developed  a  new  approach  to  longwave  computations, 
whereby  many  parameter  selections  are  automated,  permitting  reliable  use  by  less  skilled 
operators.  This  integrated  program  approach,  called  the  LongWave  Propagation  Capa¬ 
bility  (LWPC),  is  now  in  the  initial  release  process  and  is  expected  to  find  broad 
application  within  the  defense  community. 

Verification  of  the  program’s  integrity  through  a  continuing  process  of  comparison 
with  existing  and  future  VLF/LF  field  measurements  is  needed.  Such  comparative  anal¬ 
yses  guide  the  development  of  the  ionospheric  model  used  in  the  LWPC,  which  must 
be  fine-tuned  to  assure  maximum  accuracy  when  applied  to  increasingly  complex  world¬ 
wide  scenarios. 

In  an  earlier  effort  by  NOSC,  VLF/LF  data  recorded  on-board  the  merchant  ship 
GTS  Callaghan  was  reduced  and  aggregated  using  an  HP-9000  computer,  and  was  sub¬ 
sequently  compared  with  predicted  data  using  the  LWPC  model.  The  comparison  was 
limited  to  four  VLF  frequencies,  and  to  daytime-only  conditions.  In  the  present  task, 
similar  analyses  of  nighttime  data  are  conducted,  again  limited  to  the  VLF  band. 
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In  the  sections  that  follow,  this  report  summarizes  certain  adjustments  recommended 
for  the  LWPC  ionospheric  model  as  the  result  of  the  comparisons,  and  gives  some  con¬ 
clusions  about  analysis  techniques  that  could  be  used  in  future  comparisons.  The  main 
body  of  the  report  is  given  in  section  II,  where  the  objectives  and  methodology  of  the 
analysis  are  stated,  along  with  discussions  of  the  data  reduction  and  propagation  pre¬ 
diction  computer  programs  utilized  in  the  analysis.  Analysis  of  the  Callaghan  nighttime 
data  is  discussed  in  section  II-5,  and  sample  comparisons  of  predicted  and  measured  data 
are  presented  in  section  II-6.  Statistical  evaluations  of  the  data  are  shown  in  section 
II-7  along  with  summary  charts  that  provide  a  technique  for  optimized  selection  of  ion¬ 
ospheric  profiles.  Final  conclusions  and  recommendations  are  discussed  in  section  !I-8. 
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SECTION  1-2 


SUMMARY  AND  CONCLUSIONS 

This  paper  documents  a  VLF/LF  Data  Analysis  task  by  the  Naval  Ocean  Systems 
Center  to  improve  the  modeling  of  the  nighttime  ionosphere  when  making  propagation 
predictions  with  the  LWPC  computer  program.  The  task  utilizes  an  extensive  database 
of  VLF  measured  data  recorded  during  the  1985-1986  trips  of  the  merchant  ship  GTS 
Callaghan  in  the  North  Atlantic  area.  By  constraining  the  Callaghan  data  to  those  pe¬ 
riods  when  both  the  ship  and  the  distant  transmitters  were  in  time  zones  consistent  with 
all-nighttime  propagation,  and  by  eliihinating  data  from  trips  outside  the  principal  area 
of  interest,  an  aggregated  set  of  recorded  data  was  assembled  for  each  frequency  of  con¬ 
cern.  Four  frequencies  were  examined:  16.0,  19.0,  21.4  and  24.0  kHz. 

Recorded  data  sets  were  graphed  as  signal  vs.  distance  plots,  computing  distance  from 
the  transmitter  for  each  ship’s  location.  The  LWPC  program  was  then  utilized  to  com¬ 
pute  signal  vs.  distance  along  a  typical  path  in  the  same  ocean  area,  and  the  predicted 
and  recorded  data  were  compared.  By  changing  the  LWPC  parameters  p  and  h',  which 
are  exponential  relationships  describing  ionospheric  electron  density  and  collision  fre¬ 
quency  profiles  that  vary  with  height,  different  propagation  predictions  were  compared 
with  the  recorded  data  until  a  best  fit  was  obtained. 

In  all,  some  265  separate  LWPC  runs  were  compared  with  measured  data  at  four 
frequencies,  using  manual  overlays  on  a  light  table  to  find  the  best  match  between  pre¬ 
dicted  and  measured  curves.  Although  this  technique  enabled  the  selection  of  improved 
p-h'  values,  the  method  was  labor  intensive  and  required  subjective  judgments  that  be¬ 
came  difficult  at  the  higher  frequencies.  However,  best-fit  ^-h'  values  for  each 
frequency  were  selected  by  that  method,  and  a  final  confirmation  was  made  using  a  sta¬ 
tistical  evaluation  technique  wherein  average  values  of  the  means  and  standard 
deviations  of  the  difference  between  recorded  and  predicted  signal  levels  were  graphed. 

The  present  LWPC  default  values  for  nighttime  p  and  h'  are  0.30-0.80  km’^  and  87 
km,  respectively.  (P  is  entered  as  a  range  corresponding  to  the  frequency  range  of  10-60 
kHz;  LWPC  then  finds  P  for  a  specific  frequency  by  interpolation.)  The  analysis  effort 
determined  that  a  p  range  of  0.40-0.90,  and  an  h'  of  85  provided  a  better  fit  between 
recorded  and  predicted  data,  with  the  improvement  becoming  quite  significant  at  the 
lower  frequencies.  A  preliminaiy  recommendation  for  LWPC  ^-h'  profiles  is  shown 
in  figure  1. 
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For  the  higher  frequencies  of  .21.4  and  24.0  kHz,  propagation  paths  from  the  CONUS 
sources  of  these  signals  fell  above  the  70°  dip  angle  line  for  a  portion  of  the  paths  (see 
figure  3).  Since  the  LWPC  program  adjusts  the  default  )3-h'  profiles  for  polar  cap 
(PCAP)  effects  on  paths  traversing  areas  above  the  70“  line,  and  since  part  of  the  task 
was  to  determine  the  validity  of  the  PCAP  default  as  well  as  that  of  the  /3-h'  profiles, 
the  analysis  of  the  higher  frequencies  became  more  complex.  To  examine  PCAP  effects, 
duplicate  LWPC  runs  were  made.  Normal  runs  utilized  the^  default  PCAP  setting  of  70- 
74,  while  excursion  runs  set  the  PCAP  range  at  80-84,  which  had  the  effect  of  moving 
the  dip  angle  line  well  above  the  propagation  paths  of  interest,  thus  preventing  automatic 
LWPC  profile  modification.  The  results  of  these  tests  showed  better  data  fits  for  the 
excursion  settings,  indicating  that  the  default  PCAP  values  may  require  adjustment  at 
the  higher  VLF  frequencies.  Results  for  both  settings  are  shown  in  figure  1. 

It  was  concluded  that  the  present  LWPC  default  settings  for  nighttime  ionospheric 
profiles  should  be  adjusted.  Preliminary  profile  recommendations  are  to  increase  the 
P  range  by  0.10  km*^  and  to  decrease  h'  by  2  km.  The  default  PCAP  settings  of  70-74 
may  also  require  adjustment,  but  further  data  analysis  is  required  before  further  sug¬ 
gestions  may  be  made  in  this  area.  The  following  recommendations  are  made: 

•  LWPC  profile  optimization  should  be  mechanized  by  modifying  the  statis¬ 
tical  graphing  program  to  automatically  determine  the  average  of  mean  and 
standard  deviation  values  of  recorded-minus-predicted  data,  for  a  number 
of  distance  ranges.  Resulting  averages  should  be  plotted  as  a  function  of 
^  and  h'  as  a  means  of  finding  truly  optimized  profiles  for  each  frequency. 

•  The  optimization  of  LWPC  m’ghttime  profiles  should  be  extended  over  a 
much  wider  range  of  ^-h'.  This  larger  effort  would  be  made  feasible  by 
the  automatic  statistical  graphing  technique  suggested  above. 

•  The  more  extensive  LWPC  runs  should  be  facilitated  by  utilizing  the  newer, 
machine-independent  LWPC,  Version  1.0,  which  will  permit  faster  runs  and 
eliminate  the  need  for  VAX- VMS  procedures. 

•  Daytime  profile  optimization  of  the  LWPC  as  determined  in  the  Callaghan 
Report  should  be  re-examined,  using  automated  statistical  profile  compar¬ 
isons  over  a  wider  /S-h'  range. 

•  Both  daytime  and  nighttime  comparisons  should  be  extended  to  higher  fre¬ 
quencies  as  measured  data  becomes  available.  At  a  minimum,  the  LF 
channels  of  the  Callaghan  database  should  be  examined. 
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FREQUENCY  (kHz) 


^  Preliminary  Q  With  -  LWPC  Default 

Recommendation  PCAP=80,84  Profile 

With  PCAP=70.74 


Figure  1.  Preliminary  recommendation, 

revised  ranges  for  LWPC  nighttime  propagation. 
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SECTION  IM 


OBJECTIVES 

The  overall  task  objective  is  to  provide  an  in-depth  comparative  analysis  between  ex¬ 
isting  NOSC  VLF/LF  field-strength  measurements  and  signal-strength  predictions  as 
generated  by  the  LWPC  computer  program.  A  principal  goal  is  to  develop  recommen¬ 
dations  for  the  specification  of  refined  LWTC  ionospheric  ifi-h')  profiles,  so  that  future 
LWPC  predictions  agree  more  closely  with  observed  data.  In  this  task,  the  analysis  is 
constrained  to  VLF-band  frequencies  and  nighttim.e  propagation  only. 

An  interim  report*  documented  the  following  task  elements: 


•  Review  documentation  and  computer  programs  which  comprise  the 
LWPC.  Perform  LWPC  test  computations  on  contractor  computers 
and  compare  results  with  similar  calculations  supplied  by  NOSC. 

•  Review  NOSC-supplied  LWPC  predictions  as  given  in  the  1987  Cal¬ 
laghan  report,  and  compare  with  daytime  propagation  data  from  the 
GTS  Callaghan  measurements.  Develop  a  standard  methodology 
for  graphic  comparison  of  Callaghan  data  with  new  LWPC  calcula¬ 
tions  generated  within  the  task.  Compare  with  1987  Callaghan 
report  results  in  order  to  verify  analysis  methodology. 

The  present  report  documents  results  from  the  principal  task  elements,  which  are: 


•  Perform  LWPC  calculations  for  nighttime  cases  and  compare  with 
measurements  furnished  by  NOSC. 

•  Provide  comparative  analyses  of  measured  data  and  LWPC  predic¬ 
tions.  Provide  recommendations  for  refinements  in  LWPC 
ionospheric  model  to  better  align  measured  and  computed  data. 


•“Comparative  Analysis  -  Callaghan  Data”,  TCS  Report  TP-89-177,  31  May  1989. 
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SECTION  11-2 


METHODOLOGY 

The  methodology  of  the  VLF/LF  Data  Analysis  task  was  constrained  by  that  of  the 
Callaghan  Report;  every  effort  is  made  to  utilize  common  methods  and  data  reduction 
techniques  to  assure  a  continuity  of  results.  The  basic  technique  consists  of  comparing 
plotted  data  from  the  Callaghan  recorded  measurements  with  similar  plots  of  LWPC  sig¬ 
nal  predictions.  The  LWPC  ionospheric  profile  parameters  ^  and  h'  are  then  adjusted 
as  necessary  to  obtain  a  best  fit  between  the  two  curves. 

Data-reduction  programs  supplied  by  NOSC  were  modified  and  used  to  produce 
curves  of  three  primary  types: 


•  Amplitude  vs.  distance  plots  of  raw  data  from  each  trip  of  the  GTS 
Callaghan. 

•  Amplitude  vs.  distance  plots  of  aggregated  data  from  trips  grouped 
by  number,  season  and/or  time-of-day. 

•  Statistical  analysis  charts  which  permit  quantitative  evaluations  of 
recorded  vs.  predicted  data  comparisons. 

Data  comparisons  and  best-fit  adjustments  of  and  h'  were  made  by  overlaying  pre¬ 
dicted-data  plots  with  those  of  aggregated  recorded  data,  using  the  analyst’s  best 
judgment.  Statistical  analysis  charts  were  then  produced  to  quantify  the  best-fit  choices. 

Analysis  of  a  typical  new  case  (a  new  transmitter  and/or  frequency)  followed  these 
procedures: 

1.  Choose  a  transmitter,  frequency,  time-of-day  and  season. 

Produce  plots  of  aggregated  recorded  data,  selecting  a 
group  of  trips  by  number  from  summary  lists  of  ship’s 
movements.  Figure  2  illustrates  a  group  of  Callaghan 
trips,  plotted  on  a  gnomonic  projection  map. 

Individual-trip  data  plots  may  be  consulted  as  necessary 
in  making  the  selection.  Alternately,  data  from  all  50 
trips  may  be  used  in  the  aggregated  plots. 
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Figure  2.  Map  showing  Callaghan  trips  selected  for  21.4  &  24,0  kHz, 
with  70“  dip  angle  limit  shown  as  wide  line. 


2.  Choose  a  propagation  path  from  the  desired  transmitter 
that  best  approximates  the  North  Atlantic  ocean  area  tra¬ 
versed  by  the  selected  ship  movements.  Figure  3* 
illustrates  the  four  propagation  paths  utilized  in  the  pres¬ 
ent  task  effort.  The  wide  line  denotes  the  position  of  the 
70°  dip  angle  limit. 

3.  Using  the  LWPC  program,  compute  signal  vs.  distance 
data  for  the  selected  path,  transmit  frequency  and  radi¬ 
ated  power  level.  Produce  a  plot  for  comparison  with  the 
aggregated  recorded  data  plot. 

4.  Iterate  the  LWPC  calculation  as  necessary  over  a  range 
of  ionospheric  profiles.  Select  best-fit  comparisons  with 
the  recorded  data  plot.  It  may  be  necessary  to  select  sev¬ 
eral  predicted  data  plots  for  later  selection  during  the 
statistical  analysis  process. 

5.  Produce  a  statistical  analysis  of  each  best-fit  selection  as 
a  plot  of  signal  histograms  vs.  distance.  Compute  the  av¬ 
erage  of  the  mean  and  standard  deviation  values  for  each 
applicable  distance.  Plot  the  averages  as  a  function  of  ^ 
and/or  h'  for  final  best-fit  selection. 


Because  certain  of  the  Callaghan  data  trip  files  were  found  to  be  unsuitable,  it  was 
necessary  to  make  a  preliminary  selection  of  recorded  data  to  be  aggregated.  Section 
II-5  discusses  the  factors  considered  in  this  selection. 


•The  map  of  figure  3  was  made  using  four  stacked  runs  of  the  NOSC  program  WMAP.  TCS 
has  modified  WMAP  to  permit  spooling  to  a  file,  with  later  off-line  printing  on  a  laser  printer. 
The  map  of  figure  2  was  generated  in  1  min.  30  sec.  on  the  AST-386  computer,  and  was  printed 
in  1  min.  14  sec. 


10 


11 


Figure  3.  Propagation  paths  used  for  signal-strength  calculation 

as  a  function  of  distance,  with  70“  dip  angle  line  shown. 


SECTION  11-3 


DATA  REDUCTION  COMPUTER  PROGRAMS 

The  analysis  project  utilized  a  number  of  data  reduction  programs  and  supporting 
files  that  were  employed  by  NOSC  in  the  previous  Callaghan  data  analysis.  In  all,  some 
80  files  were  received  in  formats  suitable  for  use  on  Hewlett-Packard  Model  9000  com¬ 
puters.  The  project  examined  all  files  and  converted  a  number  of  the  Fortran  programs 
for  execution  on  MS-DOS  computers. 

The  following  NOSC  codes  were  converted  and  tested  for  routine  use  by  the  project: 

•  FORMIT  Forms  a  binary  data  file  for  a  specific  trip  and  frequency.  Gen¬ 

erates  an  optional  plot  of  signal  vs.  distance. 

•  PRHEAD  Reads  first  and  last  records  of  FORMIT-produced  binary  files. 

•  PLOT_DIS  Produces  aggregated  signal  data  at  a  single  frequency  of  selec¬ 

ted  trips  for  daytime,  nighttime  or  day-night  transition  periods. 
Plots  resulting  signal  data  as  a  function  of  distance.  Requires 
individual  trip  files  from  FORMIT  as  inputs. 

•  TRIPSTAT  Produces  a  statistical  analysis  of  the  difference  between  aggre¬ 

gated  recorded  data  and  predicted  data  vs.  distance.  Samples 
at  user-selected  intervals  (now  125  km)  are  analyzed,  with  re¬ 
sults  plotted  as  histograms  of  the  distribution  from  0-dB 
difference,  plus  plots  of  the  mean  of  that  distribution  and  the 
level  exceeded  by  90%  of  the  data.  Permits  user-selected 
constraints  similar  to  those  of  PLOT__DIS,  and  requires 
FORMIT  trip  files  and  predicted  signal  levels  as  inputs.. 

Several  data-reduction  program  changes  and  enhancements  were  implemented  for 
use  in  the  project.  These  changes  included: 

•  In  FORMIT,  input  was  modified  from  an  interactive  method  to  that 
of  a  constant  namelist  filename  (FORMIT.DAT). 

•  In  FORMIT,  automatic  extraction  of  ship’s  start  and  end  dates  and 
times  from  the  CNV  files  was  added. 
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•  In  PLOT_DIS  and  TRIPSTAT,  the  SUBDB  and  DELTDB  refer¬ 
ences  were  deleted. 

•  In  PLOT.DIS,  the  redundant  KFREQ  parameter  was  deleted  from 
namelist  and  the  function  was  automatically  generated  from  the 
FRQ  input.  Missing  ordinate  and  abscissa  labels  were  added. 

•  In  FORMIT,  PLOT.DIS  and  TRIPSTAT,  laser  printer  and  land¬ 
scape  switches  were  added. 

•  In  FORMIT,  PLOT.DIS  and  TRIPSTAT,  optional  output  spooling 
to  file  was  added,  permitting  later  plotting  on  a  laser  printer  in  the 
batch  mode,  with  print  enhancements  by  user  (linewidth,  line  den¬ 
sity,  size,  x-y  position). 

•  In  PLOT_DIS,  capability  for  automatic  plotting  of  predicted  data 
from  a  FASTMC.ASC  file,  atop  recorded  data  plot,  was  added. 

•  In  TRIPSTAT,  capability  for  automatic  input  of  predicted  data  from 
a  FASTMC.ASC  file  was  added. 

•  As  initially  used,  PLOT_DIS  and  TRIPSTAT  permitted  ambiguous 
and  undetected  trip  number  errors.  For  example,  with  FORMIT 
data  files  for  trips  1-20  in  the  working  directory,  a  call  for  data  from 
trips  1-10  or  trips  11-20  yielded  plots  of  trip  1-10  data,  but  were 
marked  as  called.  The  programs  have  now  been  corrected. 

•  As  received,  PLOT^DIS  and  TRIPSTAT  were  capable  of  handling 
data  from  only  20  and  19  trips,  respectively.  This  limitation  has 
been  removed. 

•  PLOT_DIS  and  TRIPSTAT  have  been  compiled  for  32-bit  opera¬ 
tion,  speeding  execution  of  each  by  more  than  two-to-one  with 
respect  to  the  original  8-bit  compiles. 


13 


The  analysis  project  also  converted  all  TRIP.NRM  files  of  recorded  data  and  the  ac¬ 
companying  CNV.DAT  files  of  ship’s  movement  from  the  52  numbered  trips,  and 
installed  the  files  on  project  computers.  Of  the  52  numbered  trips  of  the  GTS  Callaghan, 
data  from  trips  35  and  42  were  found  to  be  missing.  Therefore,  nighttime  analysis  was 
based  on  data  from  the  remaining  50  trips. 

A  listing  of  the  50  trips  of  the  GTS  Callaghan  is  given  in  table  1,  on  the  next  page. 
The  start  and  end  dates  and  times  (in  Universal  Time)  differ  slightly  from  those  listed 
in  table  i  of  the  Callaghan  report,  for  trips  1  through  22.  The  dates  and  times  of  table 
1  were  taken  directly  from  the  CNV  files  of  ship’s  movement,  and  were  utilized  in  the 
analyses  unless  recorded  data  near  the  origin  and/or  destination  were  found  to  be  in¬ 
appropriate. 
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Table  1.  GTS  Callaghan  Trips 


TRIP 

ORIGIN 

DESTINATION 

START 

END 

1 

BAYONNE,  SC 

BREMMERHAVEN 

03-29-85 

2230 

04-06-85 

1700 

2 

BREMERHAVEN.W.QR. 

CHARLESTON 

04-11-85 

2200 

04-1945 

0900 

3 

CHARLESTON.  SC 

BAYONNE 

04-20-85 

0500 

04-2145 

0700 

4 

BAYONNE 

BREMERHAVEN 

04-24-85 

0400 

054245 

1018 

5 

BREMERHAVEN 

BAYONNE 

05-05-85 

0030 

05-1345 

0830 

6 

BAYONNE 

BREMERHAVEN 

05-15-85 

0242 

05-2245 

0800 

7 

BREMERHAVEN 

BAYONNE 

05-23-85 

2100 

05-3045 

0800 

8 

BAYONNE 

BREMERHAVEN 

05-31-85 

2300 

064745 

0530 

9 

BREMERHAVEN 

BAYONNE 

06-06-85 

2000 

06-1445 

2000 

10 

BAYONNE 

BREMERHAVEN 

06-17-85 

0400 

06-2345 

0800 

11 

BREMERHAVEN 

BAYONNE 

0^25-85 

0400 

074145 

0730 

12 

BAYONNE 

CHARLESTON 

07-02-85 

0100 

07-0345 

0700 

13 

charleston 

SAVANNAH 

07-04-85 

0500 

07-0445 

0930 

14 

SAVANNAH 

ALEXANDRIA  EGYPT 

07-08-85 

0430 

07-1945 

0230 

15 

ALEXANDRIA 

BREMERHAVEN 

07-21-85 

1400 

07-2745 

2130 

16 

BREMERHAVEN 

BAYONNE 

07-29-85 

0000 

0844-85 

1715 

17 

BAYONNE 

BREMERHAVEN 

06436-85 

0100 

08-14-85 

2200 

18 

BREMERHAVEN 

ALEXANDRIA 

08-16-85 

0200 

08-22-85 

0330 

19 

ALEXANDRIA 

SAVANNAH 

06-25-85 

0700 

094645 

0712 

20 

SAVANNAH.  QA 

BREMERHAVEN 

09-1(M5 

0330 

09-1845 

1900 

21 

BREMERHAVEN 

NORFOLK 

09-20-85 

2400 

09-2845 

2000 

22 

NORFOLK.  VA 

CHARLESTON 

09-30-85 

0119 

09-30-85 

1900 

23 

NORFOLK 

CHARLESTON 

10-23-85 

1600 

10-2445 

0800 

24 

CHARLESTON 

ROHERDAM,  HOLLAND 

10-25-85 

1300 

114145 

1700 

25 

ROHERDAM 

LIVERPOOL 

11-03-85 

1445 

11-0445 

2000 

26 

LIVERPOOL  ENGLAND 

ANTWERP 

11-064)5 

1534 

114745 

1722 

27 

ANTWERP.  BELGIUM 

BREMERHAVEN 

11-0645 

2400 

11-0945 

0830 

28 

BREMERHAVEN 

CHARLESTON 

11-12-85 

0300 

11-2145 

0830 

29 

CHARLESTON 

BAYONNE 

11-22-85 

0500 

11-2345 

1100 

30 

BAYONNE 

NORFOLK 

11-2845 

0130 

11-2845 

1000 

31 

NORFOLK 

SAVANNAH 

1241-85 

2300 

12-0245 

1630 

32 

SAVANNAH 

CRISTOBAL  PANAMA 

12-0545 

14» 

12-0645 

0641 

33 

CRISTOBAL 

SAVANNAH 

12-1045 

1740 

12-1345 

0630 

34 

SAVANNAH 

PUERTO  CORTEZ.  HONDURAS 

12-1745 

0200 

12-1945 

1122 

35 

NO  DATA 

36 

CRISTOBAL 

JACKSONVILLE 

12-2345 

2214 

12-2645 

1100 

37 

JACKSONVILLE.  FL 

ZEEBRUGGE 

12-2945 

0700 

01-0646 

1646 

38 

ZEEBRUGGE,  BELGIUM 

BREMERHAVEN 

01-0846 

0654 

01-0646 

1925 

39 

BREMERHAVEN 

BAYONNE 

01-1046 

2300 

01-1946 

1300 

40 

BAYONNE 

CHARLESTON 

02-04-86 

2330 

02-0646 

1116 

41 

CHARLESTON 

ROHERDAM 

02-0846 

0400 

02-1646 

0154 

42 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-2246 

0600 

034446 

0718 

44 

BEAUMONT,  TX 

CHARLESTON 

034746 

1830 

034846 

0918 

45 

CHARLESTON 

ROHERDAM 

034946 

0600 

03-1646 

0500 

46 

ROHERDAM 

BOGEN.  NORWAY 

03-18-86 

0300 

03-2046 

1000 

47 

BOGEN 

MOREHEAD  CITY.  NC 

03-24-86 

0015 

044246 

1300 

48 

MOREHEAD  CITY 

CHARLESTON 

044546 

0015 

04-0546 

0900 

49 

CHARLESTON 

BAYONNE 

0446-86 

0500 

0447-86 

0830 

50 

BAYONNE 

ROHERDAM 

044846 

0100 

04-15-86 

0900 

51 

ROHERDAM 

BREMERHAVEN 

04-16-86 

0230 

04-1646 

1101 

52 

BREMERHAVEN 

BAYONNE 

04-18-86 

0100 

04-2546 

1700 
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SECTION  11-4 


LWPC  COMPUTER  PROGRAM 

LWPC,  Version  0,  as  written  for  the  VAX- VMS  system,  has  been  converted  for  use 
on  MS-DOS  computers  by  effort  in  this  project  and  in  associated  programs  for  the  De¬ 
fense  Nuclear  Agency.  LWPC  was  compiled  for  32-bit  operation  on  80386  machines 
using  the  Lahey  F-77L/EM-32  Fortran  compiler,  and  executes  at  speeds  comparable  to 
or  exceeding  those  of  the  DEC  Micro VAX-II  computer.  Specific  LWPC  program  mod¬ 
ules  available  for  project  use  are: 

•  PRESEG 

•  MODEFNDR 

•  SEGMWVGD 

•  MF.SW 

•  MERGE 

•  FASTMC 

•  PLTFMC 

LWPC  is  executed  on  MS-DOS  computers  using  the  VCL  (Virtual  Command  Lan¬ 
guage)  program,  which  creates  an  environment  similar  to  that  of  VAX  computers  using 
DCL  (DEC  Command  Language).  Thus,  the  various  DCL  command  files  created  by 
LWPC  in  its  automatic  sequences  may  be  used  and  interpreted  without  change.  Oper¬ 
ation  follows  descriptions  in  the  NOSC  LWPC  User’s  Guide,  as  supplemented  by  a  con- 
tractor-written  User’s  Guide  for  MS-DOS  computers. 

Use  of  the  LWPC  programs  was  validated  in  two  ways.  First,  a  reference-path  cal¬ 
culation  was  produced  using  input  data  from  the  LWPC  User’s  Guide.  A  plot  from  this 
calculation  was  determined  to  be  equivalent  to  figure  D-3  of  the  NOSC  User’s  Guide. 
Second,  the  LWPC  was  used  to  compute  signal  vs.  distance  data  for  the  four  curves  of 
VLF  predicted  data  shown  in  the  Callaghan  Report.  Computations  for  these  four  paths 
were  compared  with  those  of  the  Callaghan  report  and  documented  in  an  earlier  project 
report*. 


*“Comparative  Analysis  -  Callaghan  Data”,  TCS  Report  TP-89-177,  31  May  1989. 
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A  number  of  LWPC  program  changes  and  enhancements  were  implemented  for  use 
in  the  project.  These  included: 


•  Automatic  operation  under  emulated  VAX-VMS  environment. 

•  Filename  conventions  revised  to  accommodate  the  8  +  3  character 
filename  limitation  of  MS-DOS. 

•  Automatic  generation  of  Norton  FILEINFO.FI  files  to  annotate  re¬ 
stricted  filenames.  Initial  implementation  adds  path  bearings  and 
date/time  tags. 

•  In  PRESEG,  four-digit  power  levels  may  now  be  specified. 

•  In  PRESEG,  the  NPRINT  function  was  added  for  automatic  feed 
to  the  FASTMC  command  file. 

•  In  FASTMC  and  PLTFMC,  Amplitude/Phase  vs.  Distance  print  for¬ 
mat  (for  NPRINT  >  0)  was  changed  to  nine  columns  (instead  of 
twelve)  for  easier  printing  on  80-column  printers. 

•  In  FASTMC,  a  two-column  listing  of  distance  and  amplitude  in  an 
automatically-generated  .ASC  file  was  added,  for  use  in  inputting 
predicted  data  automatically  to  PLOT_DIS  and  to  TRIPSTAT. 

•  In  PLTFMC,  NRPTS  was  added  to  namelist  input  parameters. 

•  In  PLTFMC,  landscape  mode  was  added. 

•  In  PLTFMC,  an  output  switch  was  added  to  permit  selection  of  plot¬ 
ter  or  laser  printer  HPGL  emulation. 

•  In  PLTFMC,  output  spooling  to  file  was  added,  permitting  later 
plotting  on  a  laser  printer  in  batch  mode,  vwth  print  enhancements 
by  the  user  (linewidth,  line  density,  size,  x-y  position,  etc.). 

•  Automatic  configuration  management  and  archiving  of  all  program 
versions  by  number  was  implemented. 
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SFXTIONII-5 

MEASURED  NIGHTTIME  DATA 
II-5.1  SELECTED  TRIP  SETS 

Because  of  the  relative  variability  of  recorded  nighttime  signal  data,  with  respect  to 
the  more  uniform  signals  of  the  daytime  environment,  it  was  important  to  constrain  the 
Callaghan  measurements  as  much  as  possible  before  comparing  with  the  LWPC  predic¬ 
tions.  Normalized  data  from  each  trip  was  examined  for  insured  integrity  at  each 
frequency  of  interest,  rejecting  those  cases  where  obvious  anomolies  in  the  measure¬ 
ments  had  occurred  or  where  the  ship’s  course  was  outside  the  area  represented  by  the 
selected  propagation  path  to  be  computed  by  the  prediction  model.  The  data  reduction 
program  PLOT_DIS  was  run  for  individual  trip  data  at  each  frequency,  and  a  reference 
manual  of  plotted  data  was  assembled  as  an  aid  to  selecting  the  groups  of  valid  trips. 

Individual  plots  of  the  ship’s  course  for  each  trip  were  also  made,  using  the  Callaghan 
navigation  files.  The  NOSC  program  WMAP  was  modified  to  permit  automatic  plotting 
of  the  ship’s  course  on  a  world  map,  and  a  reference  manual  of  individual  trip  plots  was 
created.  (Figure  2  is  an  example  of  such  plots  for  a  group  of  trips.) 

A  “Primary  Trip  Set”  of  recorded  data  was  then  selected,  following  the  procedure 
diagramed  in  figure  4,  and  repeating  the  process  for  each  of  the  four  frequencies.  Trip 
data  files  containing  other  than  all-nighttime  data  were  first  eliminated.  Plots  of  the 
ship’s  course  for  each  trip  were  then  examined,  and  trips  out  of  the  North  Atlantic  area 
of  principal  interest  were  rejected,  along  with  several  short  trips  in  the  European  area 
and  between  CONUS  ports.  Some  of  the  trips  to  Central  America  were  retained  for 
the  16.0-  and  19.0-kHz  British  transmitters,  smce  such  trips  appeared  to  provide  useful 
data  in  the  6-  to  8-Mm  range.  The  Central  America  trips  were  rejected  for  the  CONUS- 
based  transmitters,  however,  since  such  locations  were  considerably  removed  from  the 
propagation  paths  of  interest.  Finally,  the  remaining  trip  data  was  examined  for  integ¬ 
rity,  rejecting  files  that  were  obviously  contaminated  by  noise. 

The  remaining  data  files  were  defined  as  the  Primary  Trip  Sets,  and  aggregated  re¬ 
corded  data  from  these  sets  were  used  for  the  comparisons.  Tables  2-5  on  the  following 
pages  record  the  final  trip  selections,  and  indicate  the  reasons  for  rejecting  individual 
data  files  using  the  reject  codes  of  figure  4. 
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REJECT  CODE  1 


REJECT  CODE  2 


REJECT  CODE  3 


REJECT  CODE  4 


Figure  4.  Methodology  diagram,  Callaghan  trip  recorded  data, 
selection  of  measured  data  for  comparison. 
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Table  2.  Primary  trip  set,  16.0  kHz  (large  type). 


REJECT  TRIP 

ORIGIN 

DESTINATION 

START 

END 

CODE 

1 

BAYONNE,  SC 

BREMMERHAVEN 

03-29^ 

2230 

044645 

1700 

2 

BREMERHAVEN.W.QR. 

CHARLESTON 

04-11.85 

2200 

04-1945 

0900 

2 

3 

CHAnUSTON,  sc 

BAYONNE 

04-30.W 

0900 

04-3149 

0700 

4 

BAYONNE 

BREMERHAVEN 

04-24-85 

0400 

054245 

1018 

S 

BREMERHAVEN 

BAYONNE 

0541545' 

0030 

05-1345 

0830 

4 

t 

■AYONNE 

BREMERHAVEN 

09-1949 

0343 

09-3349 

0900 

4 

7 

BRCMEWAVCN 

BAYONNE 

09-3349 

3100 

09-3049 

0900 

4 

1 

BAYONNE 

BREMERHAVEN 

094149 

3300 

090749 

0930 

9 

BREMERHAVEN 

BAYONNE 

0646-85 

2000 

06-1445 

2000 

10 

BAYONNE 

BREMERHAVEN 

06-17-85 

0400 

06-2345 

0800 

11 

BREMERHAVEN 

BAYONNE 

06-2545 

0400 

074145 

0730 

2 

•J 

BAYONNE 

CHARLESTON 

07-0349 

0100 

07-0349 

0700 

1 

13 

CHAnUSTON 

UVANNAH 

074449 

0900 

07-0446 

0930 

14 

SAVANNAH 

ALEXANDRIA.  EGYPT 

07-0645 

0430 

07-1945 

0230 

15 

ALEXANDRIA 

BREMERHAVEN 

07-21-85 

1400 

07-2745 

2130 

16 

BREMERHAVEN 

BAYONNE 

07-2945 

0000 

0844-65 

1715 

17 

BAYONNE 

BREMERHAVEN 

084845 

0100 

08-1445 

2200 

2 

•1 

BREMENHAVEN 

ALEXANDRA 

09-1946 

0300 

09-3349 

0330 

19 

ALEXANDRIA 

SAVANNAH 

06-2545 

0700 

094645 

0712 

20 

SAVANNAH,  QA 

BREMERHAVEN 

09-1045 

0330 

09-1645 

1900 

21 

BREMERHAVEN 

NORFOLK 

09-2045 

2400 

09-2848 

2000 

2 

33 

NOWOU^VA 

CHARUSTON 

09-3049 

0111 

09-3949 

1300 

2 

33 

NOPtfOlX 

CHARUSTON 

10-3349 

1900 

193449 

0900 

24 

CHARLESTON 

ROHEROAM,  HOLLAND 

10-2545 

1300 

114145 

1700 

2 

39 

ROTTEROAM 

UVERPOOL 

114349 

1449 

11-0449 

1000 

2 

3« 

irvERROOl,  ENCKANO 

ANTWERP 

114949 

1934 

11-0749 

1733 

2 

3» 

ANTWERP.  BEIWUM 

BREMERHAVEN 

11-0949 

3400 

114949 

0930 

26 

BREMERHAVEN 

CHARLESTON 

11-1245 

0300 

11-2145 

0830 

2 

3t 

CHARLESTON 

BAYONNE 

11-3349 

0900 

11-3349 

1100 

2 

30 

BAYONNE 

NORPOIK 

11-3949 

0130 

11-3946 

1000 

1 

31 

NORROIK 

SAVANNAH 

134149 

3300 

130349 

1930 

32 

SAVANNAH 

CRISTOBAL*.  PANAMA 

12-0545 

1430 

12-0645 

0641 

33 

CRISTOBAL 

SAVANNAH 

12-1045 

1740 

12-1345 

0630 

34 

SAVANNAH 

PUERTO  CORTEZ,  HONDURAS 

12-1745 

0200 

12-1948 

1122 

1 

33 

NO  DATA 

36 

CRISTOBAL 

JACKSONVILLE 

12-2345 

2214 

12-2645 

1100 

37 

JACKSONVILLE.  FL 

ZEEBRUGQE 

12-2945 

0700 

01-0646 

1846 

1 

3« 

ZEEBnJOQE,  BEiaiUM 

BREMERHAVEN 

014949 

0994 

014949 

1339 

39 

BREMERHAVEN 

BAYONNE 

01-1046 

2300 

01-1946 

1300 

2 

«o 

BAYONNE 

CHARLESTON 

034449 

3330 

034949 

1119 

1 

41 

CHARLESTON 

ROTTERDAM 

084949 

0400 

03-1349 

0194 

1 

«3 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-2246 

0600 

034446 

0718 

2 

u 

BEAUMONT,  TX 

CHARLESTON 

034749 

1930 

034949 

0919 

45 

CHARLESTON 

ROTTERDAM 

034946 

0800 

03-1646 

0500 

2 

*1 

ROTTERDAM 

BOOEN,  NCPWAV 

03-1949 

0300 

33-3049 

1800 

47 

BOGEN 

MOREHEADCITY.NC 

03-24-86 

0015 

044246 

1300 

2 

U 

MOflEhEAOaiY 

CHARLESTON 

044549 

0019 

3J14949 

0900 

2 

49 

CHARLESTON 

BAYONNE 

044949 

0900 

3^749 

0930 

50 

BAYONNE 

ROTTERDAM 

044646 

0100 

04-15-86 

0900 

1 

5’ 

ROTTERDAM 

BREMERHAVEN 

04-1649 

0330 

09-1646 

1101 

52 

BREMERHAVEN 

BAYONNE 

04-1846 

0100 

04-25-86 

1700 

20 


Table  3.  Primary  trip  set,  19.0  kHz  (large  type). 


REJECT  TRIP 

ORIGIN 

DESTINATION 

START 

END 

CODE 

1 

BAYONNE,  SC 

BREMMERHAVEN 

03-29^ 

2230 

04-0645 

1700 

2 

BREMERHAVEN.W.QR. 

CHARLESTON 

04-11-8S 

2200 

04-19-85 

0900 

2 

3 

CHARLESTON,  SC 

BAYONNE 

04.2MS 

0300 

34-31-33 

0700 

4 

4 

BAYONNE 

BREMERHAVEN 

04-34-tS 

0400 

03-03-33 

1013 

5 

BREMERHAVEN 

BAYONNE 

OS0545 

0030 

05-1345 

0630 

6 

BAYONNE 

BREMERHAVEN 

05-15^5 

0242 

05-2245 

0800 

7 

BREMERHAVEN 

BAYONNE 

05-23^ 

2100 

05-3045 

0800 

8 

BAYONNE 

BREMERHAVEN 

05-31-85 

2300 

064745 

0530 

9 

BREMERHAVEN 

BAYONNE 

06-0845 

2000 

06-1445 

2000 

10 

BAYONNE 

BREMERHAVEN 

06-17-85 

0400 

06-2345 

0600 

11 

BREMERHAVEN 

BAYONNE 

06-25-85 

0400 

0741-85 

0730 

2 

12 

BAYONNE 

CHARUSTON 

07-03-3S 

0100 

07-03-33 

0700 

1 

13 

CHARUSTON 

SAVANNAH 

07-044S 

0300 

07-04-33 

0330 

14 

SAVANNAH 

ALEXANDRIA,  EQY.OT 

07-08-85 

0430 

07-1945 

0230 

15 

ALEXANDRIA 

BREMERHAVEN 

07-21-85 

1400 

07-2745 

2130 

16 

BREMERHAVEN 

BAYONNE 

07-2945 

0000 

08-04-85 

1713 

17 

BAYONNE 

BREMERHAVEN 

084645 

0100 

06-1445 

2200 

2 

11 

BREMERHAVEN 

ALEXANDRIA 

03-134S 

0300 

03-33-33 

0330 

19 

ALEXANDRIA 

SAVANNAH 

08-2545 

0700 

094645 

0712 

20 

SAVANNAH,  QA 

BREMERHAVEN 

09-1045 

0330 

09-1845 

1900 

21 

BREMERHAVEN 

NORFOLK 

09-2045 

2400 

09-28-85 

2000 

2 

32 

NORFOLK,  VA 

CHARUSTON 

03-3&38 

0113 

03-3043 

1300 

2 

33 

NORFOLK 

CHARUSTON 

10-33-4S 

1300 

10-3443 

0300 

24 

CHARLESTON 

ROHERDAM,  HOLLAND 

10-2545 

1300 

114145 

1700 

2 

33 

ROTTERDAM 

UVfRROOL 

11-3348 

1443 

114443 

3000 

2 

M 

LWERROOL.  INOLANO 

ANTWERP 

It-OB-tS 

1334 

114748 

1733 

2 

37 

AN1WERR,  BSLOIUM 

BREMERHAVEN 

11-03« 

3400 

1143-33 

0190 

28 

BREMERHAVEN 

CHARLESTON 

1M245 

0300 

11-2145 

0830 

2 

33 

charleston 

BAYONNE 

n-33-IS 

0300 

11-3343 

1100 

2 

30 

BAYONNE 

NORFOLK 

11-3346 

0130 

11-3343 

1000 

1 

31 

NORFOLK 

SAVANNAH 

13-OI-4a 

3300 

1343-33 

1330 

32 

SAVANNAH 

CRISTOBAL,  PANAMA 

12-0545 

1430 

12-0645 

0641 

33 

CRISTOBAL 

SAVANNAH 

12-1045 

1740 

12-13-85 

0630 

34 

SAVANNAH 

PUERTO  CORTEZ,  HONDURAS 

12-1745 

0200 

12-1945 

1122 

1 

33 

NO  DATA 

36 

CRISTOBAL 

JACKSONVILLE 

12-2345 

2214 

12-2645 

1100 

37 

JACKSONVILLE,  FL 

ZEEBRUGQE 

12-2945 

0700 

01-0648 

1646 

1 

33 

ZEEBRUOOE,  BELOIUM 

BREMERHAVEN 

014346 

0334 

014343 

1333 

39 

BREMERHAVEN 

BAYONNE 

01-1086 

2300 

01-1948 

1300 

2 

40 

BAYONNE 

CHARUSTON 

034443 

3330 

034343 

1113 

1 

41 

CHARLESTON 

ROITEROAM 

034343 

0400 

03-1343 

0134 

1 

43 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-2246 

0600 

034446 

0718 

2 

u 

BEAUMONT,  TX 

CHARUSTON 

034743 

liSO 

034343 

0313 

45 

CHARLESTON 

ROHERDAM 

034946 

0800 

03-1646 

050<' 

2 

43 

ROTTERDAM 

BOOEN.  NORV/AY 

03-1343 

0300 

03-3043 

1000 

47 

BOGEN 

MCREHEADCITY.NC 

03-24-86 

0015 

0442-86 

1300 

2 

44 

MOREHEAOOTY 

CHARUSTON 

044343 

0013 

0A4543 

0900 

2 

40 

CHARLESTON 

BAYONNE 

04-03-33 

0300 

0447-M 

0130 

50 

BAYONNE 

ROHERDAM 

0448-96 

0100 

04-15-86 

0900 

1 

SI 

ROTTERDAM 

BREMERHAVEN 

04-1346 

0330 

04-13-36 

1101 

52 

BREMERHAVEN 

BAYONNE 

04-1846 

0100 

04-2546 

1700 

21 


Table  4.  Primary  trip  set,  21.4  kHz  (large  type) 


REJECT  TRIP 

ORIGIN 

DESTINATION 

START 

END 

CODE 

1 

BAYONNE,  SC 

BREMMERHAVEN 

03-29^5 

2230 

04436-85 

1700 

2 

BREMERHAVEN,  W.  QR. 

CHARLESTON 

04-11-85 

2200 

04-19-85 

0900 

2 

3 

CHAHLESTON,  SC 

BAYONNE 

0<-3049 

0900 

04-2149 

0700 

4 

BAYONNE 

BREMERHAVEN 

04-24-M 

0400 

05432-85 

1018 

9 

BREMERHAVEN 

BAYONNE 

05-09^5 

0030 

05-13-85 

0630 

6 

BAYONNE 

BREMERHAVEN 

05-15-85 

0242 

05-22-85 

0800 

7 

BREMERHAVEN 

BAYONNE 

05-23-89 

2100 

05-30-85 

0800 

8 

BAYONNE 

BREMERHAVEN 

05-31-85 

2300 

06-07-85 

0530 

9 

BREMERHAVEN 

BAYONNE 

08-08-85 

2000 

06-14-85 

2000 

10 

BAYONNE 

BREMERHAVEN 

08-17-85 

0400 

06-23-85 

0800 

11 

BREMERHAVEN 

BAYONNE 

08-25-85 

0400 

07-01-85 

0730 

2 

’2 

BAYONNE 

CHARUSTON 

07.02-I9 

01M 

07-03-95 

07W 

1 

13 

CHARLESTON 

SAVANNAH 

07-04.99 

0900 

07-04-99 

09X 

14 

SAVANNAH 

ALEXANDRIA.  EGYPT 

07-08-89 

0430 

07-19-85 

0230 

15 

ALEXANDRIA 

BREMERHAVEN 

07-21-89 

1400 

07-27-85 

2130 

18 

BREMERHAVEN 

BAYONNE 

07-29-85 

0000 

08-04-85 

1715 

17 

BAYONNE 

BREMERHAVEN 

08438-85 

0100 

08-14-85 

2200 

2 

ti 

BREMERHAVEN 

ALEXANDRIA 

09-1949 

0300 

a9-33-U 

03X 

19 

ALEXANDRIA 

SAVANNAH 

08-25-85 

0700 

09436-85 

0712 

20 

SAVANNAH,  OA 

BREMERHAVEN 

09-10-85 

0330 

09-18-85 

1900 

21 

BREMERHAVEN 

NORFOLK 

09-20-89 

2400 

09-28-85 

2000 

2 

23 

NORFOLK.  VA 

CHARLESTON 

09-30-99 

0119 

09-X49 

1900 

2 

33 

NORFOLK 

CHARUSTON 

10-3349 

1900 

10-34-99 

0900 

24 

CHARLESTON 

ROTTERDAM,  HOLLAND 

10-25-85 

1300 

11431-85 

1700 

2 

29 

ROTTERDAM 

UVERROOL 

11-0349 

1449 

11-0449 

3000 

2 

2« 

UVERROOU  INOLANO 

ANTWERP 

11-0949 

1934 

11-4749 

1733 

2 

27 

ANTWERP.  BCiaiUM 

BREMERHAVEN 

11-0949 

3400 

11-09-99 

09X 

28 

BREMERHAVEN 

CHARLESTON 

11-12-89 

O'JOO 

11-21-85 

0830 

2 

31 

charleston 

BAYONNE 

11-3349 

0900 

11-33-99 

1IM 

2 

X 

BAYONNE 

NORFOUL 

11-39-99 

01X 

11.M49 

1000 

1 

31 

NORFOLK. 

SAVANNAH 

13-0149 

3300 

13-03-99 

19X 

3 

33 

SAVANNAH 

CRISTOBAU  PANAMA 

13-09-99 

14X 

13-0949 

0941 

3 

33 

CRISTOBAL 

SAVANNAH 

13-10-99 

1740 

13-13-99 

09X 

3 

34 

SAVANNAH 

PUERTO  CORTEZ.  ‘•CNOURa 

1H749 

0300 

1M9-99 

1133 

1 

X 

NO  DATA 

3 

at 

CRISTOBAL 

JACKSONVKU 

12-2349 

3314 

13-39-99 

1100 

37 

JACKSONVILLE,  FL 

ZEEBRUGGE 

12-29-85 

0700 

01438-86 

1846 

1 

M 

ZEEBRUOQE.  BELGIUM 

BREMERHAVEN 

01-09-99 

0994 

01-0949 

1939 

39 

BREMERHAVEN 

BAYONNE 

01-10-88 

2300 

01-19-88 

1300 

2 

40 

BAYONNE 

CHARUSTON 

034449 

33X 

03-0949 

1119 

1 

4t 

CHARUSTON 

ROTTERDAM 

03-0949 

0400 

03-1949 

0194 

1 

42 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-22-88 

0600 

03434-88 

0718 

2 

44 

BEAUMONT,  tX 

CHARUSTON 

034749 

tsx 

034949 

0919 

45 

CHARLESTON 

nOHERDAM 

034)9-88 

0800 

03-16-86 

0500 

2 

49 

ROnEROAM 

BOOEN,  NORWAY 

03-19-99 

C3CO 

03-20-99 

1000 

47 

BOGEN 

MOREHEAO  CITY,  NC 

03-24-86 

0015 

044)2-88 

1300 

2 

49 

MOREHEAO  CnV 

CHARUSTON 

044949 

MIS 

04-09-99 

0900 

2 

49 

CHARUSTON 

BAYCNNE 

044949 

09X 

0447-99 

09X 

50 

BAYONNE 

ROHERDAM 

04-08-86 

O’CO 

04-15-86 

0900 

1 

3t 

ROTTERDAM 

BREMERHAVEN 

04-1949 

0233 

04-19-99 

1101 

52 

BREMERHAVEN 

BAYONNE 

04-18-88 

0100 

04-25-86 

1700 

22 


Table  5.  Primary  trip  set,  24.0  kHz  (large  type) 


REJECT  TRIP 

ORIGIN 

DESTINATION 

START 

END 

CODE 

1 

BAYONNE,  SC 

BREMMERHAVEN 

03-29^ 

2230 

044645 

1700 

2 

BREMERHAVEN,W.3R. 

CHARLESTON 

04.11^ 

2200 

04-1945 

0900 

2 

3 

CHAnLESTON,  sc 

BAYONNE 

04-3048 

0800 

04-3148 

07K 

4 

BAYONNE 

BREMERHAVEN 

04-2445 

0400 

054245 

1018 

5 

BREMERHAVEN 

BAYONNE 

054545 

0030 

0813-85 

0830 

6 

BAYONNE 

BREMERHAVEN 

05-1545 

0242 

082245 

0800 

7 

BREMERHAVEN 

BAYONNE 

05-2345 

2100 

083045 

0800 

8 

BAYONNE 

BREMERHAVEN 

054145 

2300 

064745. 

0530 

9 

BREMERHAVEN 

BAYONNE 

064845 

2000 

081445 

2000 

10 

BAYONNE 

BREMERHAVEN 

06-1745 

0400 

082345 

0600 

11 

BREMERHAVEN 

BAYONNE 

06-2545 

0400 

074145 

0730 

2 

12 

•AYONNE 

CHARLESTON 

OT-03-38 

OIK 

07-0348 

07K 

1 

13 

CHARLESTON 

SAVANNAH 

074448 

08H 

07-0440 

08K 

14 

SAVANNAH 

ALEXANDRIA,  EGYPT 

07-0845 

0430 

07-1945 

0230 

15 

ALEXANDRIA 

BREMERHAVEN 

07-2145 

1400 

07-2745 

2130 

IS 

BREMERHAVEN 

BAYONNE 

07-2945 

0000 

084445 

1715 

17 

BAYONNE 

BREMERHAVEN 

084845 

0100 

081445 

2200 

2 

11 

SREMERHAVIN 

ALSXANDRM 

03-1348 

etK 

08-3348 

03X 

10 

ALEXANDRIA 

SAVANNAH 

082545 

0700 

094645 

0712 

20 

SAVANNAH,  QA 

BREMERHAVEN 

081045 

0330 

081845 

1900 

21 

BREMERHAVEN 

NORF01J< 

082045 

2400 

082845 

2000 

2 

33 

NORFOLK,  VA 

CHARLESTON 

034048 

0113 

08-3048 

13K 

2 

33 

NORFOLK 

CHARUSTON 

10-2348 

13M 

10-3448 

03K 

24 

CHARLESTON 

ROHERDAM,  HOLLAND 

182545 

1300 

114145 

1700 

2 

2S 

ROTTEROAM 

UVERROOL 

11-0348 

1448 

11-0448 

SOK 

2 

3« 

UVERROOl.  ENQIANO 

ANTWERR 

11-03-X 

1834 

114748 

1733 

2 

27 

ANTWERR.  •ELOIUM 

BREMERHAVEN 

11-0348 

S4K 

114848 

08X 

28 

BREMERHAVEN 

CHARLESTON 

11-1245 

0300 

11-2145 

0830 

2 

33 

CHARUSTON 

BAYONNE 

11-2348 

08K 

11-3040 

11K 

2 

X 

SAYONNE 

NORFOLK 

11-1348 

01X 

11-3348 

10K 

1 

31 

NORFOLK 

SAVANNAH 

134148 

33K 

134348 

13X 

3 

32 

SAVANNAH 

CRISTOBAL.  RANAMA 

134848 

14X 

134348 

0841 

3 

X 

CR1STOEAL 

SAVANNAH 

1MMS 

1740 

1M048 

08X 

3 

34 

SAVANNAH 

RUERTO  CORTEZ.  HONDURAS 

13-1748 

08K 

1M840 

1133 

1 

X 

NO  DATA 

3 

X 

CRISTOtAL 

JACKSONVEU 

1^3348 

3314 

1^3848 

11H 

37 

JACKSONVILLE.  FL 

ZEEBRUGGE 

12-2945 

0700 

01-0646 

1846 

1 

X 

ZEESRUOOE.  KLOUM 

BREMERHAVEN 

014338 

0884 

014343 

18K 

39 

BREMERHAVEN 

BAYONNE 

01-1046 

2300 

01-1886 

1300 

2 

40 

•AYONNE 

CHARUSTON 

034448 

23X 

034848 

1113 

1 

41 

CHARLErON 

ROTTEROAM 

034348 

04K 

0M848 

0184 

1 

43 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-2246 

0600 

034446 

0718 

2 

44 

BEAUMONT,  TX 

CHARUrON 

034748 

13X 

034348 

M13 

45 

CHARLESTON 

ROTTEROAM 

034946 

0600 

081646 

0500 

2 

4« 

ROTTERDAM 

BOQEN,  NORWAY 

03-1348 

03K 

03-3040 

10K 

47 

BOGEN 

MOREHEADCITY.NC 

03-24-86 

0015 

044246 

1300 

2 

43 

MOREHEADCITY 

CHARUSTON 

04-0S-K 

K13 

04-0848 

03K 

2 

43 

CHARLESTON 

BAYONNE 

044343 

OSK 

044748 

03X 

50 

BAYONNE 

ROHEROAM 

044886 

0100 

04-1546 

0900 

1 

51 

ROTTERDAM 

BREMERHAVEN 

04-ie-M 

02X 

04-1848 

1101 

52 

BREMERHAVEN 

BAYONNE 

04-1846 

0100 

04-2546 

1700 

II-5.2  SEASON.^L  EFFECTS 


An  attempt  was  made  to  separate  the  recorded  data  into  groups,  sorted  by  season. 
Table  6  shows  the  50  Callaghan  trips  divided  into  the  three-month  seasons  of  spring, 
summer,  autumn  and  winter,  with  the  four  periods  defined  by  the  summer  solstice  on 
June  21  and  the  winter  solstice  on  December  22.  Similarly,  six-month  “warm”  and 
“cold”  seasons  were  defined  around  the  same  dates. 

It  was  decided,  however,  that  the  Callaghan  database  was  too  small  to  support  mean¬ 
ingful  separation  by  season.  Figures  5-8  illustrate  the  recorded  data  from  the  Rugby 
transmitter  broadcasting  at  16  kHz,  separated  into  spring,  summer,  autumn  and  winter 
groupings,  respectively,  as  an  example  case.  Overlays  of  these  plots  revealed  similar  pat¬ 
terns  of  recorded  signal  with  distance -there  were  no  pronounced  seasonal  effects  that 
would  justify  the  use  of  different  ionospheric  models  in  the  LWPC. 

Similar  conclusions  were  reached  in  comparing  the  six-month  seasonal  data.  Again, 
the  great  variability  of  nighttime  data  due  to  undefined  changes  in  the  propagation  chan¬ 
nel  seemed  to  obscure  any  effects  that  were  purely  seasonal. 
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WARM  I  COLD  SEASON  i  WARM  SEASON 


Table  6.  Callaghan  trips  divided  by  season. 


TRIP 

ORIGIN 

DESTINATION 

START 

END 

1 

BAYONNE.  SC 

BREMMERHAVEN 

03-29-85 

2230 

04-06-85 

1700 

2 

BREMERHAVEN,  W.  GR. 

CHARLESTON 

04-11.85 

2200 

04-19-85 

0900 

3 

CHARLESTON,  SC 

BAYONNE 

04-20-85 

0500 

04-2145 

0700 

4 

BAYONNE 

BREMERHAVEN 

04-24-85 

0400 

0542-85 

1018 

5 

BREMERHAVEN 

BAYONNE 

05-05-85 

0030 

05-1345 

0830 

6 

BAYONNE 

BREMERHAVEN 

05-15-85 

0242 

05-22-85 

0800 

7 

BREMERHAVEN 

BAYONNE 

05-23-85 

2100 

05-30-85 

0800 

8 

BAYONNE 

BREMERHAVEN 

05-31-85 

2300 

06-0745 

0530 

9 

BREMERHAVEN 

BAYONNE 

06-08-85 

2000 

06-14-85 

2000 

10 

BAYONNE 

BREMERHAVEN 

06-17-85 

0400 

06-23-85 

0800 

*1 

BREMERHAVEN 

BAYONNE 

06-25-85 

0400 

07-01-85 

0730 

BAYONNE 

CHARLESTON 

07-02-85 

0100 

07-03-85 

0700 

13 

CHARLESTON 

SAVANNAH 

07-04-85 

0500 

07-04-85 

0930 

14 

SAVANNAH 

ALEXANDRIA  EGYPT 

07-08-85 

0430 

07-19-85 

0230 

15 

ALEXANDRIA 

BREMERHAVEN 

07-21-85 

1400 

07-27-85 

2130 

'S 

BREMERHAVEN 

BAYONNE 

07-29-85 

0000 

08-04-85 

1715 

17 

BAYONNE 

BREMERHAVEN 

08-08-85 

0100 

08-14-85 

2200 

18 

BREMERHAVEN 

ALEXANDRIA 

08-16-85 

0200 

08-22-85 

0330 

19 

ALEXANDRIA 

SAVANNAH 

08-25-85 

0700 

0946-85 

0712 

20 

SAVANNAH.  QA 

BREMERHAVEN 

09-10-85 

0330 

09-18-85 

1900 

21 

BREMERHAVEN 

NORFOLK 

09-20-85 

2400 

09-28-85 

2000 

22 

NORFOLK.  VA 

CHARLESTON 

09-30-85 

0119 

09-3045 

1900 

23 

NORFOLK 

CHARLESTON 

10-23-85 

1600 

10-24-85 

0800 

24 

CHARLESTON 

ROHERDAM.  HOLLAND 

10-2545 

1300 

1141-85 

1700 

25 

ROTTERDAM 

LIVERPOOL 

11-0345 

1445 

1144-85 

2000 

25 

LIVERPOOL.  ENGLAND 

ANTWERP 

11-0645 

1534 

11-07.85 

1722 

27 

ANTWERP.  BELGIUM 

BREMERHAVEN 

11-06-85 

2400 

11-0945 

0830 

28 

BREMERHAVEN 

CHARLESTON 

11-12-85 

0300 

11-21-85 

0830 

29 

CHARLESTON 

BAYONNE 

11-22-85 

0500 

11-23-85 

1100 

30 

BAYONNE 

NORFOLK 

11-28-85 

0130 

11-28-85 

1000 

31 

NORFOLK 

SAVANNAH 

12-0145 

2300 

1242-85 

1830 

32 

SAVANNAH 

CRISTOBAL  PANAMA 

12-05-85 

1430 

12-08-85 

0641 

33 

CRISTOBAL 

SAVANNAH 

12-1045 

1740 

12-13-85 

0630 

34 

SAVANNAH 

PUERTO  CORTEZ.  HONDURAS 

12-1745 

0200 

12-19-35 

1122 

35 

NO  DATA 

36 

CRISTOBAL 

JACKSONVILLE 

12-2345 

2214 

12-2645 

1100 

37 

JACKSONVILLE.  FL 

ZEEBRUQGE 

12-29-85 

0700 

01-06-86 

1846 

38 

ZEEBRUGQE.  BELGIUM 

BREMERHAVEN 

01-0846 

0654 

01-0846 

1925 

39 

BREMERHAVEN 

BAYONNE 

01-10-86 

2300 

01-19-86 

1300 

40 

BAYONNE 

CHARLESTON 

02-04-86 

2330 

02-06-86 

1118 

-1 

CHARLESTON 

ROHEROAM 

02-08-86 

0400 

02-16-86 

0154 

42 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-22-86 

0600 

0344-86 

0718 

44 

BEAUMONT,  TX 

CHARLESTON 

03-0746 

1830 

03-08-86 

0918 

45 

CHARLESTON 

ROTTERDAM 

0349-86 

0800 

03-16-36 

0500 

45 

ROTTERDAM 

SOGEN.  NORWAY 

03-18-86 

0300 

03-20-86 

1000 

47 

EOGEN 

MOREHEAD  CITY.  NC 

03-24-86 

0015 

0442-36 

1300 

48 

MOREHEAD  CITY 

CHARLESTON 

04-05-86 

0015 

04-05-86 

0900 

49 

CHARLESTON 

BAYONNE 

04-06-86 

0500 

0447-86 

0630 

50 

BAYONNE 

ROHERDAM 

04-0846 

0100 

04-15-86 

0900 

51 

ROHERDAM 

BREMERHAVEN 

04-16-86 

0230 

04-16-36 

1101 

52 

BREMERHAVEN 

BAYONNE 

04-18-86 

0100 

04-25-86 

1700 
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TCS  PLOTJDIS  09/26/89  16:27 

Freq:  16.000  kfreq:  160 

night  (05-May-85)  to  (04-Aug-85)  deldbl;  8.  deldb2:  12. 
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TCS  PLOTJDIS  09/26/89  16:32 

Freq:  16:000  kfreq:  160 

NIGHT  (OB-Aug-85)  to  (Ol-Nov-85)  deldbl:  8.  deldb2:  12. 
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Figure  7.  Typical  recorded  data  for  Autumn  season,  Rugby  at  16.0  kHz. 


II-5.3  SOLAR  EVENT  ANALYSIS 

It  was  recognized  that  unusual  solar  activity  could  have  occurred  during  the  Callag¬ 
han  trips,  resulting  in  ionospheric  and/or  geomagnetic  field  disturbances  and  consequent 
anomolous  signal  strength  recordings.  To  investigate  this  possibility,  solar  activity  re¬ 
cords  were  obtained  from  the  Space  Environment  Service  Center  of  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA)  in  Boulder,  for  the  periods  of  the 
Callaghan  trips. 

Three  solar  phenomena  are  known  to  affect  VLF  propagation: 

1.  Ionization  due  to  X-rays  from  solar  flares  (a  daytime-only  effect). 

2.  Ionization  due  to  Solar  Proton  Events  (SPEs),  which  are  similar  to  Polar-Cap 
.‘\bsorption  (PCA)  events  and  mostly  effect  propagation  North  of  60“  N.  latitude. 

3.  Magnetic  field  disturbances  (geomagnetic  storms)  caused  by  the  flow  of  solar 
particles  into  the  earth’s  magnetosphere,  where  consequent  disturbances  are 
sometimes  severe  enough  to  affect  the  ionosphere.  Severe  geomagnetic  storms 
often  force  the  auroral  zone  South  of  its  normal  60“  N.  position. 

Little  solar  activity  was  anticipated  during  the  1985-1986  period  of  the  Callaghan 
trips,  because  the  solar  cycle  was  near  a  minimum  between  the  usual  eleven-year  peaks 
of  activity.  Only  four  SPEs  occurred  during  the  trip  period,  and  only  three  of  these  co¬ 
incided  with  actual  trip  times.  Table  7  lists  the  start  and  peak  times  (again  in  UT)  for 
applicable  SPEs.  Because  SPE  activity  was  short  in  each  case,  no  correlation  against 
measured  data  effects  was  found. 

Only  five  trips  occurred  during  times  of  significant  geomagnetic  storm  activity.  The 
five  trips,  along  with  relative  levels  of  storm  activity  are  listed  in  table  8.  Storm  levels 
varied  from  “minor”  to  “severe”  in  intensity,  where  NOAA  24-hour  geomagnetic  “A 
Index”  levels  are  defined  as: 

•  Minor  Storm:  A  Index  from  30  to  50 

•  Major  Storm:  A  Index  from  50  to  100 

* 

•  Severe  Storm:  A  Index  >  100 

Only  one  severe  storm  was  experienced  during  the  1985-1986  period,  but  it  occurred 
during  Trip  41,  which  had  been  rejected  from  the  primary  trip  sets  at  all  four  frequencies. 
The  other  storm  levels  were  judged  to  be  too  low  and  infrequent  to  correlate  with  VLF 
measured  data. 
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Table  7.  Solar  proton  events  occurring  during  Callaghan  trips. 


TRIP 

ORIGIN 

DESTINATION 

TRIP  START 

TRIP  END 

SPE  START 

SPE  PEAK 

4 

BAYONNE 

BREMERHAVEN 

04-24^5 

0400 

0502-85 

1018 

04-2505 

1430 

04-2605 

0600 

14 

SAVANNAH 

ALEXANDRIA,  EGYPT 

07.06-85 

0430 

07-19-85 

0230 

07-0905 

0235 

07-0905 

0325 

40 

BAYONNE 

CHARLESTON 

02-04.86 

2330 

02-06-86 

1118 

02-0606 

0825 

02-07.86 

1730 

- 

- 

- 

- 

- 

- 

- 

02-1406 

1155 

02-15-86 

0400 

Table  8.  Callaghan  trips  during  geomagnetic  storms  (selected  trips  in  large  type). 


STORM 

LEVEL 

TRIP 

ORIGIN 

DESTINATION 

START 

END 

MAJOR 

4 

BAYONNE 

BREMERHAVEN 

04-2405 

0400 

050205 

1018 

MINOR 

14 

SAVANNAH 

ALEXANDRIA.  EGYPT 

07-0605 

0430 

07-1905 

0230 

MIN./MAJ. 

17 

BAYONNE 

BREMERHAVEN 

08-0805 

0100 

06-1405 

2200 

MIN./MAJ. 

37 

JACKSONVILLE,  a 

ZEEBRUGGE 

12-2905 

0700 

01-0606 

1846 

SEVERE 

41 

CHAnjESTON 

nOTTSnOAM 

02-0M4 

0400 

02-1444 

0154 
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II-5.4  PCAP  EVALUATIONS  . 


Figure  2,  on  page  9,  illustrates  typical  Callaghan  trips  with  respect  to  the  position 
of  the  70°  dip  angle  line  (the  wide  line).  Those  trips  which  pass  North  of  Scotland  are 
seen  to  cross  an  extensive  ocean  area  North  of  the  line.  The  70°  dip  angle  line  is  also 
shown  on  figure  3  with  respect  to  the  four  signal  propagation  paths  selected  for  predicted 
data.  Transmissions  from  the  two  CONUS  sources  at  21.4  and  24.0  kHz  pass  North  of 
the  70°  line  in  the  initial  portion  of  their  paths.  The  transmissions  from  the  British  trans¬ 
mitters,  however,  at  16.0  and  19.0  kHz,  pass  below  the  70°  line. 

The  nighttime  ionosphere  is  known  to  be  more  depressed  in  the  polar  regions,  be¬ 
coming  more  day-like  at  the  higher  latitudes.  The  LWPC  program  recognizes  this  effect 
by  including  a  polar-cap  (PCAP)  transition  that  modifies  the  /5-h'  model  across  a  user- 
selectable  boundary.  The  default  PCAP  setting  is  70-74,  corresponding  to  a  transition 
area  from  70°  to  74°  in  dip  angle,  ^-h'  profiles  are  automatically  adjusted  for  propa¬ 
gation  paths  passing  through  the  PCAP  transition.  Table  9  is  an  example  listing  of 
segments  along  the  60°  propagation  path  from  the  Annapolis  transmitter,  where  Rho  is 
distance  along  the  path.  As  the  path  moves  through  the  area  North  of  the  70°  dip  angle 
line,  and  the  codip  angle  falls  below  20°  (codip  =  90°  -  dip),  a  corresponding  change  is 
seen  in  ^-h'.  Table  10  lists  the  Callaghan  trips  that  traversed  areas  below  the  70°  dip 
angle  line,  and  identifies  those  selected  for  the  primary  trip  sets  at  21.4  and  24.0  kHz. 

The  Callaghan  data  is  thus  a  mixed  bag  from  the  standpoint  of  polar-cap  effects,  and 
it  proved  difficult  to  isolate  PCAP-free  data.  Instead,  the  project  ran  LWPC  in  the  de¬ 
fault  setting  (PCAP  =  70-74)  for  all  frequencies,  then  added  excursions  at  21.4  and  24.0 
kHz  where  PCAP  was  set  at  80-84,  thus  moving  the  transition  area  well  above  the  prop¬ 
agation  path  areas.  Results  of  the  PCAP  excursions  are  presented  in  the  next  sections. 


Table  9.  Path  segmentation  for  Annapolis,  path  bearing  60°. 


RHO 

AZIM 

CODIP 

MAGFLD 

SIGMA 

EPSR 

BETA 

HPRIME 

0.000 

068.0 

19.4 

0.542 

4e+00 

81.0 

0.39 

84.8 

0.060 

069.1 

19.2 

0.542 

3e-03 

15.0 

0.39 

84.8 

0.260 

072.8 

18.7 

0.543 

4e-00 

81.0 

0.39 

84.8 

0.800 

083.1 

18.0 

0.539 

4e-00 

81.0 

0.38 

82.7 

1.180 

090.0 

18.0 

0.533 

O 

o 

•1 

a> 

81.0 

0.39 

84.8 

2.480 

106.5 

20.0 

0.502 

4e-00 

81.0 

0.41 

87.0 

5.920 

120.0 

30.0 

0.435 

<D 

O 

ro 

15.0 

0.41 

87.0 

6.240 

120.8 

31.3 

0.431 

CM 

o 

Q> 

15.0 

0.41 

87.0 

6.460 

121.4 

32.3 

0.428 

4e-r00 

81.0 

0.41 

87.0 

7.060 

123.1 

35.2 

0.421 

4e+00 

81.0 

0.41 

87.0 
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Table  10.  Callaghan  Trips  below  the  70*  dip  angle  line  (large  type). 
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SECTION  11-6 

PREDICTED  VS.  MEASURED  DATA  COMPARISONS 
II-6.1  RUGBY,  16.0  kHz 

Primary  results  of  the  analysis  are  presented  in  this  and  the  following  sections.  As 
was  described  in  the  methodology  discussion  of  section  II-2,  PLOT.DIS  amplitude  vs. 
distance  plots  of  aggregated  data  from  the  Callaghan  recordings  are  compared  with  over¬ 
lays  of  predicted  data  graphed  to  the  same  scale.  Figures  9-38  are  plots  of  such  data 
arranged  as  recorded  and  predicted  pairs,  with  the  measurements  and  predictions  shown 
as  a  function  of  distance  from  the  transmitter  along  great  circle  paths.  All  data  is  for 
nighttime,  extracting  trip  data  satisfying  the  criteria  that  both  the  transmitter  and  re¬ 
ceiver  (ship)  locations  have  solar  zenith  angles  greater  than  90“. 

Graphs  for  the  Rugby  transmitter  at  16.0  kHz  are  given  in  figures  9-13,  which  are 
illustrative  of  the  many  comparisons  made  during  the  analysis.  Because  of  the  great  vol¬ 
ume  of  such  graphs  generated  by  the  project,  only  representative  samples  are  presented 
in  this  report.  A  common  format  is  used  for  all  frequency  sections,  where  the  first  graph 
compares  recorded  data  with  the  LWPC  default-setting  prediction.  A  sequence  of  four 
additional  comparisons  then  follows,  illustrating  a  different  selection  spread  over  four 
values  of  h'.  The  discussions  of  section  II-7.7  e.xplain  how  summary  graphs  were  used 
to  verify  ^-h'  profile  recommendations  for  the  overall  frequency  range. 

The  16-kHz  examples  of  this  section  indicate  how  the  recorded  vs.  predicted  data 
fit  was  improved  by  increasing  and  decreasing  h'  with  respect  to  the  LWPC  default 
settings.  The  comparisons  concentrate  on  the  distance  ranges  between  zero  and  five 
megameters  in  searching  for  best  data  fits.  It  is  possible  to  improve  the  fit  in  the  ex¬ 
treme  far  field  (5-8  Mm)  by  increasing  substantially,  but  the  close-in  curves  no  longer 
match  as  well  and  the  selection  becomes  suspect.  Constraint  of  the  data-fit  selection 
to  distances  <  5  Mm  is  justified  because,  beyond  that  range,  the  measured  data  is  less 
reliable  as  recorded  levels  approach  those  of  atmospheric  and  shipboard  noise. 

In  each  of  the  figures,  applicable  trip  numbers  are  shown  below  the  lower  edge  of 
the  graph.  The  predicted  data  file  noted  above  the  graph  gives  a  two-character  run  code, 
then  the  letters  “as”  denoting  the  LWPC  ASCII  output  file,  followed  by  the  )3-h'  dves- 
ignaiioa.  For  example,  the  notation  of  figure  9  shows  that  "Run  fc”  was  for  a  range 
of  0.30  to  0.80  (corresponding  to  a  10-  to  60-kHz  frequency  range)  and  was  for  an  h' 
of  87  km. 
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TCS  PLOTJDIS  09/26/89  15:09  Predicted  data  file:  fcas3080.B70 

Freq:  16.000  kfreq:  160 

NIGHT  {29-Mar-B5)  to  (25-Apr-86)  deldbl:  8.  deldbP:  12. 


(uj/Ani  0AoqB  gp)  apnqnduiv 


Distance  (Mm) 


TCS  PLOTJDIS  09/26/89  19:  10  Predicted  data  file:  eeas4O9O.05O 

Freq:  16.000  kfreq:  160 

.NIGHT  (29-Mar-85)  to  (25-Apr-86)  deldbl;  8.  de3db2:  12. 


(ui/Ani  aAoqe  gp)  apnqiiduiv 


Distance  (Mm) 


II-6.2  ANTHORNE,  19.0  kHz 

Results  for  the  19.0-kHz  transmitter  are  shown  in  figures  14-18.  As  was  the  case 
for  16  kHz,  the  default  profile  shows  the  longer  “cycles”  between  fades  when  compared 
with  measured  data,  indicating  that  a  change  to  a  lower  h'  is  in  order.  This  is  confirmed 
in  figures  15-18,  where  ^  was  adjusted,  also. 

/3  is  increased  to  the  0.40-0.90  range  in  this  instance.  With  interpolation  by  the 
LWPC,  a  ^  of  0.49  was  selected  by  the  computer  for  this  frequency. 
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CS  PLOTJDIS  09/26/89  15:27  Predicted  data  file:  fcas3080.B70 

req:  19.000  kfreq:  190 

NIGHT  (29-Mar-85)  to  (25-Apr-86)  deldbl:  8.  deldb2;  12. 


(ui/An^  aAoqe  gp)  apn^ncliuv 


Anthorne 


CS  PL0TJ3IS  09/26/89  15:31  Predicted  data  file:  geas4090.B40 

req:  19.000  kfreq:  190 

NIGHT  (29-Mar-85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12. 


(lu/Ant  8Aoq9  gp)  apn:induiv 


LWPC  prediction  for  ^  =  0.49,  h  =84,  PCAP= 70-74, 
with  recorded  data  for  Anthorne  at  80  kW,  19.0  kHz. 


TCS  PL0TJ3IS  09/26/89  15:36  Predicted  data  file:  eeas4090.850 

Freq:  19.000  kfreq:  190 

•  NIGHT  (29-Mar-a5)  to  (25-Apr-8B)  deldbl;  8.  deldb2:  12 


(lu/AOT  8Aoqe  gp)  apnqndiuv 


Distance  (Mm) 


CS  PLOTJDIS  09/26/89  15:40  Predicted  data  file:  heas4090.860 

req:  19.000  kfreq:  190 

NIGHT  (29-Mar-a5)  to  (25-Apr-86)  deldbl:  8.  deldb2;  12. 


LWPC  prediction  for  ^=0.49,  h'  =  86,  PCAP  =  70-74, 
with  recorded  data  for  Anthorne  at  80  kW,  19.0  kHz, 


II-6.3  ANNAPOUS,  21.4  kHz  ^ 

As  frequency  increases,  wavelength  shortens  and  the  earth-ionosphere  waveguide  can 
support  many  more  propagation  modes.  This  is  confirmed  by  the  added  complexity  of 
both  measured  and  predicted  data  at  21.4  kHz,  when  compared  with  the  similar  curves 
of  the  lower  frequencies. 

The  default  LWPC  profile  shown  in  figure  19  reflects  the  addition  of  higher-order 
modes,  but  lacks  the  deep  fading  patterns  suggested  by  the  first  two  nulls  of  the  mea¬ 
sured  data.  These  nulls  are  restored  in  the  ^3=0.51  profiles  of  figures  20-23,  but  the 
selection  of  the  best  h'  has  become  more  difficult  and  subjective  at  this  higher  frequency. 
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CS  PLOTJJIS  09/26/89  16:03  Predicted  data  file:  geas4O9O.06O 

req:  21.400  kfreq:  214 

NIGHT  (30-Mar-85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12. 


LWPC  prediction  for/?=0.51,  h'=86,  PCAP  =  70-74, 
with  recorded  data  for  Annapolis  at  250  kW,  21.4  kHz. 


CS  PLOTJDIS  09/26/89  19:23  Predicted  data  file:  ieas4090.870 

req:  21.400  kfreq:  214 

NIGHT  {30-Mar-85)  to  (25-App-B6)  deldbl;  8.  deldb2:  12. 


{ui/AHT  8Aoqe  gp)  apnqndiuv 


Distance  (Mm) 


II-6.4  CUTLER,  24.0  kHz 

The  default  LWPC  prediction  of  figure  24  shows  even  more  high-order  mode  effects 
than  in  the  21.4-kHz  case.  Again,  the  nulls  seem  less  pronounced  than  in  the  measured 
data. 

Selection  of  a  higher  ^  range,  as  in  figures  25-28,  results  in  a  partial  improvement, 
but  the  match  of  fading  positions  is  disappointing.  This  sequence  suggests  the  need  for 
automation  of  the  selection  process  so  that  many  more  profiles  may  be  examined  within 
practical  time  constraints.  Mechanization  of  the  comparisons  is  discussed  in  section  II- 
7.7. 

The  predicted  curves  shown  thus  far  in  the  figures  utilized  the  LWPC  default  polar- 
cap  settings  (see  discussion  in  section  II-5.4).  As  a  further  point  of  interest,  PCAP 
adjustments  were  eliminated  in  the  excursions  of  the  following  sections. 
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PL0TJ3IS  09/26/89  19;  51  Predicted  data  file:  icas30B0.B70 

q;  24.000  kfreq:  240 

GHT  (29-Mar-85)  to  (25-Apr-86)  deldbl;  8.  deldb2:  12. 


(lu/An;  aAoqe  gp)  apntjndmv 


LWPC  prediction  for/?  =  0.44,  h'  =  87  (default),  PCAP  =  70-74, 
with  recorded  data  for  Cutler  at  1000  kW,  24.0  kHz. 


II-6.5  EXCURSIONS,  ANNAPOLIS,  21.4  kHz,  WITH  PCAP  =  80-84 

The  graphs  of  figures  29-33  repeat  the  profile  selections  of  section  II-6.3,  except  for 
the  elimination  of  the  LWPC’s  polar-area  adjustment  of/3-h'  by  setting  PCAP  =  80-84. 
Comparison  of  the  two  sets  of  graphs  shows  the  introduction  of  even  more  high-order 
modes,  and  the  best-fit  selection  is  still  a  difficult  choice. 

One  trend  is  evident:  For  this  selection  of  deletion  of  the  PCAP  adjustment  per¬ 
mits  a  better  match  to  measured  data  at  a  lower  h'.  This  group  of  graphs  favors  the 
curve  fit  of  h'  =  85,  whereas  in  the  PCAP  =  70-74  cases  a  higher  h'  is  preferred. 
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CS  PLOTJDIS  09/26/89  19:28  Predicted  data  file:  deas3080.870 

req:  21.400  kfreq:  214 

NIGHT  (30-Mar-85)  to  (25-Apr-8B)  deldbl:  8.  deldb2:  12. 


LWPC  prediction  for  fi  =  0Al,  h'  =  87,  (default),  PCAP 
with  recorded  data  for  Annapolis  at  250  kW,  21.4  kHz. 


Predicte'd  data  file;  laas4090.B40 


LWPC  prediction  for /3  =  0.51,  h  =84, 

with  recorded  data  for  Annapolis  at  250  kW,  21.4  kHz. 


CS  PL0TJ3IS  09/26/89  19:41  Predicted  data  file:  lbas4090.B60 

req:  21.400  kfreq:  214 

NIGHT  (3O-Mar-05)  to  (25-Apr-8B)  deldbl:  8.  deldb2:  12. 


LWPC  prediction  for /3  =  0.51.  h  =86,  PCAP  =  80-84, 
with  recorded  data  for  Annapolis  at  250  kW,  21.4  kHz. 


CS  PL0TJ3IS  09/26/89  19:46  Predicted  data  file:  lcas4090.870 

req:  21.400  kfreq:  214 

NIGHT  OO-Mar-BS)  to  (25-Apr-8B)  deldbl:  8.  deldb2:  12. 


Distance  (Mm) 


II-6.6  EXCURSIONS,  CUTLER,  24.0  kHz,  WITH  PCAP  =  80-84 

Similar  tradeoffs  were  noted  in  the  24.0-kHz  cases  of  figures  34-38.  The  selection 
of  a  higher  ^  range  and  the  trend  to  lower  h'  values  resulted  in  better  curve  fits  than 
in  the  default  PCAP  graphs. 

No  final  conclusions  could  be  reached  about  the  PCAP  adjustment  in  these  limited 
tests;  data  fit  seemed  to  suggest  a  tradeoff  of  PCAP  range  vs.  h'.  NOSC  has  suggested 
that  more  extensive  tests  of  PCAP  settings,  taken  over  a  much  wider  y3-h'  range,  should 
be  conducted,  and  have  recommended  at  least  three  ranges  for  the  PCAP  transition  anal¬ 
ysis:  80-84,  70-74  and  66-70. 


S  PLOT_DIS  09/26/89  21:44  Predicted  data  file:  lbas4090.860 

eq:  24.000  kfreq:  240 

IGHT  (29-Mar-Q5)  to  (25-Apr-8B)  deldbl;  8.  deldb2:  12. 


LWPC  prediction  for  ^  =  0.54,  h'  =  86,  PCAP  =  80-84, 
with  recorded  data  for  Cutler  at  1000  kW,  24.0  kHz. 


SECTION  11-7 


STATISTICAL  ANALYSES 


II-7.1  RUGBY,  16.0  kHz 

While  the  graphs  of  the  previous  section  are  useful  for  quick  comparisons  of  the  total 
data  set  by  the  human  eye,  the  NOSC  program  TRIPSTAT  provides  a  better  technique 
for  quantitative  evaluation.  In  section  II-7,  statistical  graphs  from  TRIPSTAT  runs  are 
presented  for  each  case  of  the.  earlier  sequence. 

The  statistical  charts  of  recorded-minus-predicted  signal  strength  in  dB  at  16.0  kHz 
are  given  in  figures  39-43.  Each  chart  presents  separate  histograms  at  125-km  intervals 
along  the  path,  evaluating  the  distribution  of  data  points  recorded  within  a  range  of  ±60 
km  from  each  125-km  interval.  The  total  number  of  recorded  points  within  each  ±60- 
km  range  are  graphed  in  the  lower  chart  of  each  sheet. 

Each  histogram  is  divided  into  1-dB  bins,  the  height  of  which  varies  in  proportion 
to  the  number  of  data  points  in  each  1-dB  interval.  In  addition,  mean  (diamond)  and 
90%  exceedance  level  (x)  values  are  displayed  below  each  histogram.  An  overall  mea¬ 
sure  of  goodness  is  the  degree  with  which  the  histogram  and  its  mean  are  centered  on 
the  0-dB  axis.  A  0-dB  mean  with  all  data  points  centered  on  the  axis  would  represent 
a  perfect  match  between  recorded  and  predicted  signal  levels  at  that  distance. 

The  statistical  graph  for  the  LWPC  default  profile  at  16  kHz  is  shown  in  figure  39, 
where  the  cyclic  nature  of  the  recorded-minus-predicted  data  reflects  the  mismatch  be¬ 
tween  the  fading  cycles  displayed  in  figure  9,  page  35.  The  distribution  of  means  is 
improved  somewhat  by  the  series  of  charts  in  figures  40-43,  but  these  quantitative  anal¬ 
yses  confirm  that  the  nighttime  ionosphere  is  quite  variable  with  respect  to  the  much 
closer  fits  of  daytime  data  in  the  earlier  Callaghan  report.  This  quantitative  evaluation 
nevertheless  permits  a  better  evaluation  of  one  profile  against  the  other,  as  is  seen  by 
simply  scanning  the  distribution  of  mean  values. 

The  statistical  graphs  were  the  main  means  for  optimizing  the  selection  of  best-fit 
profiles  in  this  task,  considering  the  limited  number  of  profiles  examined.  A  method 
for  automating  this  selection  over  a  greater  number  of  choices  is  discussed  in  section 
II-7.7. 
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TCS  TRIPSTAT  09/17/09 

Predicted  data  file:  eeas4O9O.05O 
Freq:  16.000  Months:  01  to 


Figure  41.  Recorded  -  predicted  data,  Rugby  a 
with  LWPC  prediction  for  p  =  0.46,  li 


II-7.2  ANTHORNE,  19.0  kHz 

The  statistical  comparison  between  the  LWPC  default  and  the  optimized  ^-h'  values 
is  more  favorable  in  the  19.0*kHz  examples  of  figures  44-48.  The  distribution  of  means 
in  the  default  case  is  scattered  over  a  wide  range,  while  the  means  of  the  85-km  example 
of  figure  46  are  relatively  well  behaved,  with  most  falling  within  a  ±4-dB  range  out  to 
4.5-Mm  distance. 

The  distribution  of  90%  exceedance  values  is  seen  as  a  less  favorable  indicator  of 
fit  than  is  that  of  the  means  in  figure  46,  where  the  existence  of  spurious  data  points, 
particularly  around  fades  as  is  typical  of  the  nighttime  ionosphere,  produces  large  values 
for  the  90%  exceedance  points. 
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II-7.3  ANNAPOUS,  21.4  kHz 

The  LWPC  default  profile  provides  a  much  better  fit  in  the  21.4-kHz  case,  as  is  seen 
by  the  charts  of  figure  49.  The  charts  of  figures  50-53,  for  the  ^  =  0.40-0.90  profile  are 
worse  than  those  of  figure  49  at  the  shorter  waveguide  heights.  The  best  fit  for  this 
profile  sequence  is  found  for  h'  =  87  km,  in  figure  53,  where  all  but  one  mean  value  fall 
within  a  ±3-dB  range,  and  this  fit  extends  past  the  5-Mm  distance. 
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TCS  TRIPSTAT  09/17/89 

Predicted  data  file:  keas4090.840 
Free:  21.400  Months:  01  to 


Figure  50.  Recorded  —  predicted  data,  Annapolis 
with  LWPC  prediction  for  ^  =  0.51,  h'  = 


II-7.4  CUTLER,  24.0  kHz 


The  statistical  evaluations  for  24.0  kHz  show  that  the  ^-h'  =  0.40-0.90  profile  offers 
an  improvement  over  the  LWPC  default  profile,  as  is  shown  in  figures  54-58.  The  dis¬ 
tribution  of  mean  values  is  again  best  at  h'  =  87  km,  but  fewer  points  are  within  the 
±3-dB  range. 
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TCS  TRIPSTAT  09/17/89  19:  22 

Predicted  data  file:  feas4090.850 

Freq:  24.000  Months:  01  to  12  NIGHT  Deldis:  60  Xmax:  500 
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Figure  56.  Recorded  —  predicted  data,  Cutler  at  1000  kW,  24.0  kHz, 
with  LWPC  prediction  for/?  =  0.54,  h'  =  85,  PCAP  =  70-74. 


TCS  TRIPSTAT  09/17/09  19:  00 

Predicted  data  file:  ieas4090.870 

Freq;  24.000  Months:  01  to  12  NIGHT  Deldis:  60  Xmax:  500 
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Figure  58.  Recorded  —  predicted  data.  Cutler  at  1000  kW,  24.0  kHz, 
with  LWPC  prediction  for  y3= 0.54,  h'  =  87,  PCAP= 70-74. 


II-7.5  EXCURSIONS,  ANNAPOUS,  21.4  kHz,  WITH  PCAP  =  80-84 

Figures  59-63  examine  the  statistical  characteristics  of  the  21.4-kHz  path  with  the 
PCAP  adjustment  of  LWPC  disabled.  As  was  seen  in  the  signal  vs.  distance  plots,  re¬ 
moval  of  the  PCAP  compensation  forces  the  best-fit  case  to  a  lower  reflection  height. 
The  best  case  for  these  examples  is  found  at  h'  =  85  km. 

In  the  85-km  charts  of  figure  61,  most  mean  values  to  beyond  5-Mm  distance  are 
within  =2  dB  of  the  zero  line.  This  agreement  between  measured  and  predicted  data 
is  indeed  remarkable,  given  the  variable  nature  of  the  nighttime  ionosphere,  and  ap¬ 
proaches  the  agreement  achieved  for  the  daytime  cases  of  the  Callaghan  report. 
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TCS  TRIPSTAT  09/19/89  09:56 

Predicted  data  file:  deas3090.B70 

Freq:  21.400  Months:  01  to  12  NIGHT  Deldis:  60  Xmax:  500 


(ap)  paTotpadd  -  papdooafcj  STunoo 
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Figure  59.  Recorded  —  predicted  data,  Annapolis  at  250  kW,  21.4  kHz, 

with  default  LWPC  prediction  for /3  =  0.41,  h'=87,  PCAP  =  80-84. 


II-7.6  EXCURSIONS,  CUTLER,  24.0  kHz,  WITH  PCAP  =  80-84 

The  excursions  of  figures  64-68,  for  24.0  kHz,  reflect  the  variability  of  the  signal  vs. 
distance  plots  seen  earlier.  Although  the  fit  isn’t  as  good  as  in  the  21.4-kHz  examples, 
the  best  case  is  still  found  at  84  to  85  km.  This  suggests  that  the  85-km  profile  be  se¬ 
lected  as  a  standard  improvement  for  all  nighttime  frequencies,  since  it  would  then  be 
a  common  choice.  (A  preliminary  recommendation  for  nighttime  profiles  from  the  com¬ 
parisons  of  this  task  was  presented  in  figure  1  on  page  5.)  However,  such  a  decision 
ignores  the  effect  of  LWPC’s  PCAP  adjustment,  which  certainly  must  remain  in  place 
at  some  setting  to  compensate  for  the  lower  ionospheres  at  higher  latimdes. 

These  findings  underscore  the  need  for  additional  analyses  of  the  Callaghan  data  over 
a  wider  range  of  profiles,  and  for  the  investigation  of  more  wide-ranging  PCAP  selec¬ 
tions.  Tne  manner  in  which  the  statistical  comparisons  could  be  automated  is  discussed 
in  next  section. 
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II-7.7  SUMMARY  OF  STATISTICAL  COMPARISONS 

The  previous  sections  showed  that  TRIPSTAT  statistical  charts  provide  a  means  for 
quantitative  evaluation  of  recorded-minus-predicted  data  that  is  superior  to  the  analyst’s 
subjective  judgment.  Still,  the  comparison  of  one  statistical  chart  against  another  re¬ 
quires  a  similar  visual  appraisal. 

As  a  way  of  further  mechanizing  the  comparisons,  an  average  of  all  the  histogram 
mean  values  was  determined  for  each  chart  over  given  distance  ranges.  Then,  a  single 
number  was  derived  for  each  profile  and  distance  range,  and  these  numbers  were 
graphed  over  a  range  of  profile  choices.  Graphs  of  the  average  of  mean  values  appear 
in  this  section,  and  represent  summaries  of  the  final  profiles  selected  for  comparison. 

Similarly,  standard  deviations  of  each  histogram  were  computed  from  the  mean  and 
90%  exceedance  levels  and  were  averaged  for  each  chart  over  a  particular  distance 
range,  assuming  a  normal  distribution  of  data.  Again,  a  single  number  corresponding 
to  average  standard  dev^<»'.ion  was  derived  for  each  profile  and  distance,  ard  these  num¬ 
bers  were  graphed  for  a  range  of  profile  choices. 

Figures  69-74  present  graphs  of  averages  of  the  mean  and  standard  deviation  values 
vs.  h'  for  the  selected  range  of  0.40-0.90.  Figures  69-72  summarize  the  four  analysis 
frequencies,  using  the  default  PCAP  settings  of  70-74.  Summaries  of  the  21.4-  and  24.0- 
kHz  excursions  with  PCAP  *  80-84  are  given  in  figures  73-74.  In  each  graph,  averages 
for  five  distance  ranges  are  given,  and  a  point  representing  the  single  average  for  the 
LWPC  default  profile  at  a  0-5  Mm  range  is  also  shown. 

Choosing  a  profile  where  the  summary  graph  average  of  means  is  at  a  minimum  is 
the  final  step  of  the  analysis  process.  For  example,  in  figure  70  for  19-kHz,  the  0-5  Mm 
curve  shows  the  lowest  average  mean  at  an  h'  just  below  85  km. 

A  scan  of  the  summary  curves  confirms  that  h'  *  85  km  is  a  good  choice  for  the 
=  0.40-0.90  profile  range,  providing  that  PCAP  =  80-84  at  the  higher  frequencies.  Ide¬ 
ally,  three-dimensional  plots  with  both  and  h'  treated  parametrically  over  a  much  wider 
profile  range  should  be  produced,  and  PCAP  adjustment  effects  should  be  investigated. 

The  summary  graphs  show  how  the  technique  for  finding  the  average  of  means  and 
standard  deviations  can  aid  in  the  optimization  process.  The  summaries  of  this  task  were 
derived  manually  from  the  TRff  STAT  charts.  Future  tasks  should  facilitate  the  process 
by  adding  program  routines  to  automatically  produce  listings  and  graphs  of  such  averages 
from  a  full  range  of  yS-h'  profiles. 


106 


Avarag#  of  Std.  Dav.  Valuas  (dB)  Avaraga  of  Mean  Values  (dB) 


83.5  84  84.5  85  85.5  86  86.5  87  87.5 

h  ’  (km) 

Distonca 

Range  (Mm):  ••  0-1  —  0-2  -0-3  —  0-4  —  0-5 

e  For  Beta  0.3-0.8,  h‘  87  (LWPC  Default)  at  0-5  Mm. 


Figure  69.  Average  of  mean  and  standard  deviation  values  vs.  h', 
recorded  -  predicted  data  for  Rugby,  65  kW,  16.0  kHz, 
for  p  range  0.40-0.90,  PCAP= 70-74. 
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Average  of  Sfd.  Dev.  Values  (dB)  Average  of  Mean  Values  (dB) 


h  ’  (km) 

Distance 

Range  (Mm):  ••  0-1  -- 0-2  -0-3  -0-4  -  0-5 

e  For  Beta  0.3-0.8,  h’  87  (LWPC  Default)  at  0-5  Mm. 


Figure  70.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Anthorne,  80  kW,  19.0  kHz, 
for  ^  range  0.40*0.90,  PCAP  =  70-74. 
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Average  of  Sid.  Dev.  Values  (dB)  Average  of  Mean  Values  (dB) 


h  •  (km) 

Distance 

Range  (Mm):  ••  0-1  --  0-2  -0-3  —  0-4.  ——  0-5 

©•  For  Beta  0.3-0. 8,  h*  87  (LWPC  Default)  at  0-5  Mm. 


Figure  71.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Annapolis,  250  kW,  21.4  kHz, 
for^  range  0.40-0.90,  PC4P  =  70-74. 
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Average  of  Std.  Dev-  Values  (dB)  Average  of  Mean  Values  (dB) 


h  *  (km) 

Distance 

Range  (Mm):  ••  0-1  -- 0-2  -0-3  —  0-A  -  0-5 

0  For  Beta  0.3-0.8,  h’  87  (LWPC  Default)  at  0-5  Mm. 


Figure  72.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Cutler,  1000  kW,  24.0  kHz, 
for  ^  range  0.40-0.90,  PCAP  =  70-74. 


Average  of  Std.  Dev.  Values  (dB)  Average  of  Mean  Values  (dB) 


h  '  (km) 

Distance 

Range  (Mm):  ••  0-1  --  0-2  —  0-3  —  0-4  —  0-5 

6-  For  Beta  0.3-0.8,  h’  87  (LWPC  Default)  at  0-5  Mm. 


Figure  73.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Annapolis,  250  kW,  21.4  kHz, 
for  p  range  0.40*0.90,  PCAP  =  80-84. 
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Values  (dB)  Average  of  Mean  Values  (dB 


h  •  (km) 

Distance 

Range  (Mm):  0-1  —  0—2  —  0-3  —  0-4  •—  0-5 

Q  For  Beta  0.3-0.8,  h’  87  (LWPC  Default)  at  0-5  Mm. 


Figure  74.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Cutler,  1000  kW,  24.0  kHz, 
for  ^  range  0.40-0.90,  PCAP  =  80-84. 


SECTION  8 


CONCLUSIONS 

All  objectives  of  Part  1  of  the  present  task  have  been  met,  and  the  following  con¬ 
clusions  have  been  made. 

•  The  analysis  methodology  of  the  Callaghan  report  has  been  repro¬ 
duced,  and  comparable  output  graphics  and  data  reduction 
techniques  have  been  reduced  to  practice.  The  NOSC  LongWave 
Prediction  Capability  (LWPC)  family  of  programs  has  been  installed 
on  task  computers  in  an  emulated  VAX-VMS  environment  and  us¬ 
age  has  been  validated  against  NOSC-supplied  test  cases. 

•  The  daytime  analyses  of  the  Callaghan  report  have  been  repeated 
with  nearly  identical  results,  as  a  further  measure  of  methodology 
validation.  Example  predictions  of  the  same  type  have  been  pro¬ 
duced  for  the  nighttime  environment,  and  initial  comparisons  with 
recorded  data  at  21,4  kHz  show  good  agreement. 

•  Statistical  summary  graphs  should  be  expanded  to  also  include  dis¬ 
tance  ranges  of  1-5,  2-5,  3-5  and  4-5  Mm. 

•  It  is  evident  that  separate  graphs  of  recorded  data  by  trip  will  be 
helpful  in  eliminating  questionable  data  and  in  optimizing  a  pre¬ 
ferred  set  of  aggregate  data  for  each  frequency.  Such  optimized  sets 
of  recorded  data  should  be  routinely  assembled  before  profile  mod¬ 
ifications  are  attempted. 

•  Because  the  presence  of  significant  high-order  mode  interaction 
tends  to  modify  the  dominant  mode  propagation  signature,  and  be¬ 
cause  this  effect  varies  with  frequency,  the  optimization  of  nighttime 
profiles  is  expected  to  be  complex.  Profile  modifications  for  the 
nighttime  VLF  environment -the  main  objective  of  this  task- 
should  therefore  be  recommended  only  after  all  available 
frequencies  are  considered. 

•  The  Callaghan  report  examined  four  VLF  frequencies:  16.0,  19.0, 

21.4  and  24.0  kHz.  Data  from  the  other  VLF  channels  of  the  re¬ 
corded  database  (16.4,  24.8  and  28.5  kHz)  should  also  be  examined 
within  the  constraints  of  time  and  funding. 

•  If  possible,  future  analyses  should  separate  recorded  data  into  areas 
of  latitude  below  and  above  the  codip  =  20**  point  where  the  present 
LWPC  design  forces  a  profile  change. 
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